High Resolution Imaging of Single Atoms in a Quantum Gas 
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Our knowledge on ultracold quantum gases is 
strongly influenced by our ability to probe these 
objects PQ [2l |3l |4] • In situ imaging combined with 
single atom sensitivity is an especially appealing 
scenario as it can provide direct information on 
the structure and the correlations of such sys- 
tems [5j |6l [7]' For a precise characterization a 
high spatial resolution is mandatory. In partic- 
ular, the perspective to study quantum gases in 
optical lattices [51 [9l HOj makes a resolution well 
below one micrometer highly desirable. Here, we 
report on a novel microscopy technique which is 
based on scanning electron microscopy and allows 
for the detection of single atoms inside a quan- 
tum gas with a spatial resolution of better than 
150 nm. Imaging a Bose-Einstein condensate in 
a one-dimensional optical lattice with 600 nm pe- 
riod we demonstrate single site addressability in a 
sub-/im optical lattice. The technique offers excit- 
ing possibilities for the preparation, manipulation 
and analysis of quantum gases. 

Ultracold atoms can be visualized by various tech- 
niques. Absorption imaging [IT] is the workhorse in most 
experiments and is typically applied in time of flight in 
order to increase the cloud size and reduce the optical 
density. While phase contrast imaging [TJ [12] is well 
suited for trapped quantum gases, fluorescence imaging 
Q~3j HU El EI El E] is especially attractive as it allows 
for single atom detection with almost 100% efficiency. It 
has been applied to isolated thermal atoms at low den- 
sities but has not yet been extended to single atom de- 
tection in quantum gases. The best achievable resolution 
of these optical techniques is ultimately limited by half 
the wavelength of the used light field - in practice, the 
best reported resolution is about 1 jam M]. Direct parti- 
cle detection of metastable atoms in time of flight [2j E] 
and outcoupling of single atoms from a condensate with a 
radio frequency field [3] are alternative techniques which 
have been developed. However, they either cannot be ap- 
plied to trapped samples [2] E] or are restricted to one 
spatial dimension [3 . Whereas each of those techniques 
has its specific advantages and applications, a versatile in 
situ detection of single atoms in a quantum gas is lack- 
ing. Moreover, a spatial resolution of below 1 /am which 
opens the intriguing perspective to resolve single sites in 
a sub-//m optical lattice has not yet been achieved. 
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FIG. 1: Working principle. The atomic ensemble is prepared 
in an optical dipole trap. An electron beam with variable 
beam current and diameter is scanned across the cloud. Elec- 
tron impact ionization produces ions which are guided with 
an ion optical system towards a channeltron detector. The 
ion signal together with the scan pattern is used to compile 
the image. 



In our experiment we have transferred the principles 
of scanning electron microscopy to the detection of ultra- 
cold atoms (Fig. 1). A focused electron beam with 6keV 
electron energy, a full width half maximum (FWHM) 
diameter of 100-150 nm and a current of 10-20 nA is 
scanned across a Bose-Einstein condensate of rubidium 
atoms which is prepared in an optical dipole trap [20 . 
The atoms are ionized by electron impact ionization, ex- 
tracted with an electrostatic field and subsequently de- 
tected by an ion detector. The small diameter of the 
electron beam ensures a high spatial resolution, whereas 
the ion detection provides single atom sensitivity. The 
total ionization cross section at 6 keV electron energy for 
rubidium is a{ on = 3.5 x 10 _17 cm 2 [21] and represents 
40% of all scattering events [22] [23] . Elastic and inelastic 
electron-atom collisions constitute the remaining events 
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and lead to atom loss with no detectable signal. As the 
cross section is eight orders of magnitude smaller com- 
pared to the absorption cross section of a resonant pho- 
ton, the atomic cloud is optically thin for the electron 
beam. For typical parameters, only 1 out of 500,000 in- 
cident electrons undergoes a collision. 




FIG. 2: Electron microscope image of a trapped Bose-Einstein 
condensate. In a we show an image of a 87 Rb condensate. 
The image has 400 x 150 pixels with a pixel size of 300 nm 
x 300 nm. Each pixel was illuminated for 2 fis with the elec- 
tron beam (140 nm FWHM beam diameter). Every dot cor- 
responds to a detected atom. In total, 350 ions were col- 
lected during the exposure. The condensate contains about 
10 5 atoms. The sum over 300 images is presented in b. Each 
image was taken in a separate experimental run. 

When exposed to the electron beam, the probability 
for the detection of an atom at a position {x, y} is given 
by (see Methods section) 

P{x,y) = ^o-ion x At x 7? det x N J dz\(j)(x,y,z)\ 2 . (1) 

Here, / is the electron beam current, e is the electron 
charge, At is the pixel dwell time of the electron beam, 
?7det is the detector efficiency, J dz | </>(#, z)\ 2 denotes 
the column density of the atom's wave function along the 
propagation direction of the electron beam (z-direction) 
and N is the number of atoms in the single particle state. 
In Fig. 2a we show a scanning electron microscope image 
of a Bose-Einstein condensate. For our experimental pa- 
rameters, a fraction of 350 atoms is detected (the total 
number of atoms in the condensate is about 100,000). In 
a Bose-Einstein condensate all atoms occupy the same 
quantum state and the many-body wave function i/j sep- 
arates into the product of N identical single particle wave 
functions ip(xi, ...,xn) = YliLi •> with N being the 

number of atoms in the condensate. Therefore, the inter- 
pretation of the image involves quantum-mechanical con- 
cepts: As the single particle wave function <j> extends over 
the whole atomic cloud, the spatially resolved detection 
of an atom must be understood as a projective measure- 
ment in position space. As a consequence, the retrieved 



image is intrinsically probabilistic. This is in contrast 
to almost all microscopy images showing the distribution 
of individual atoms as in these cases the location of the 
atoms is already fixed prior to their detection. Another 
important aspect is related to the Heisenberg uncertainty 
principle. During the detection process, the atom is cou- 
pled to a probe (in our case an electron beam) and energy 
as well as momentum can be exchanged between them. 
Consequently, the localization of an atom within a range 
Ax enforces a momentum spread of Ap > h/(2Ax). If 
Ax is smaller than the extension of the wave function </>, 
substantial momentum transfer is unavoidable and the 
detected atom is no longer part of the condensate, regard- 
less of the specific experimental realization. Hence, the 
ionization of the atoms in our scheme does not constitute 
a serious limitation or drawback. It is even advantageous 
because it helps rapidly extracting the reaction products 
from the remaining system, keeping possible perturba- 
tions small. 

Whether the image in Fig. 2a is indeed a probabilis- 
tic selection of the full atomic distribution according to 
equation ([!]), can be checked by summing over many im- 
ages (Fig. 2b) and comparing them to a theoretical den- 
sity profile. The profile is derived from the so-called semi- 
ideal model [24 , 25 , 26 , which describes a bimodal distri- 
bution at finite temperature. While the condensate part 
is obtained from a numerical solution of the 3D-Gross- 
Pitaevskii equation, the thermal component is modelled 
as a non-interacting gas in an effective potential, taking 
into account the repulsion of the thermal atoms by the 
condensed atoms (see Methods section). The comparison 
with our data (Fig. 3a, b) exhibits very good agreement 
over the whole extension of the cloud including the wings 
of thermal atoms. This gives not only indirect evidence 
of the repulsion between the condensate fraction and the 
thermal component in the trap (Fig. 3c), but also con- 
firms that the image shown in Fig. 2a displays a proba- 
bilistic selection of atoms. 

Comparing the condensate with and without exposure 
to the electron beam (absorption images in the inset of 
Fig. 3a) we do not find any significant difference apart 
from a reduction in atom number by about 7%. These 
losses are composed of two contributions: primary elec- 
tron atom collisions and secondary collisions of the pri- 
mary reaction products. We find that every scattered 
atom or produced ion kicks off on average one more atom. 
In all these collisions the energy transfer is much larger 
than the depth of the optical potential and all scattered 
particles can escape from the trap. Essentially no energy 
is deposited in the cloud as we observe an additional 
heating of merely 5 nK after exposure to the electron 
beam. Thus, the perturbation caused by the detection 
process is very small. If not, the scanning speed could 
be made larger than the speed of sound in the conden- 



FIG. 3: Analysis of the Bose-Einstein condensate. In a we 
present the axial distribution of the condensate shown in 
Fig. 2b, integrated in the radial direction. The experimental 
data (blue columns) are compared to a bimodal distribution 
(red line) calculated in the semi-ideal model for a total atom 
number of N — 115, 000 and a temperature of T — 80 nK, cor- 
responding to a condensate fraction of 80%. The inset shows 
absorption images of the condensate after 15 ms time of flight 
with and without exposure to the electron beam. The num- 
ber of atoms after exposure is reduced by 7% in average. The 
distribution in the radial direction is shown in b. In c we 
have plotted the radial density of the thermal component in 
the trap centre as calculated from the model. The minimum 
is due to the repulsion from the condensate fraction. 



sate providing an effectively unperturbed cloud during 
the whole imaging sequence. According to equation 0, 
high imaging speed is associated with a reduced signal 
and a convenient setting of the imaging parameters has 
to be chosen for each application. Most detected ions 
are singly charged (80%) but we also find higher charged 
states of up to Rb T+ resulting from inner shell ionization. 
Only 1 out of 50 detected events is due to background gas 
ionization or dark counts which results in a high signal to 
noise ratio as evidenced by Fig. 2. Taking into account 
a detector efficiency of 30% the total efficiency for our 
detection scheme is currently limited to 12%. It can be 
increased by a more efficient ion detector and additional 
photoionization of inelastically scattered atoms. We es- 
timate that a total detection efficiency of more than 50% 
could be feasible. 

In order to characterize the resolution of our imag- 
ing technique we have loaded the condensate in a one- 
dimensional optical lattice with 600 nm lattice period. A 
sequence of electron microscope images with increasing 
resolution is shown in Fig. 4a-c. The periodic structure 
of the potential is clearly resolved with high contrast. 
As a scanning probe technique is used for the image for- 



FIG. 4: Images of a Bose-Einstein condensate loaded in a 
one-dimensional optical lattice. The lattice which is created 
by two laser beams (A = 850 nm) intersecting under 90° has a 
period of / = 600 nm. Each image is the sum of 50 individual 
images. The pixel size is 200 nm x 200 nm (a), 75 nm x 75 nm 
(b), and 25 nm x 25 nm (c). The lattice depth was 20 recoil 
energies E r (E r — 7r 2 h 2 /(2ml 2 ), with m being the rubidium 
mass) and the FWHM diameter of the electron beam was 
95 nm. 

mation, the addressability of individual lattice sites is 
demonstrated as well. The atomic density in each lattice 
site is radially symmetric with a diameter of 6 /im and 
a thickness of 300 nm and documents the large depth of 
focus of the electron optical imaging system. 

One of the most intriguing properties of a Bose- 
Einstein condensate is its macroscopic phase coherence. 
In a periodic potential the phase coherence can be easily 
verified by interference experiments. An absorption im- 
age of the condensate after a ballistic expansion of 15 ms 
is shown in the inset of Fig. 5a. The image was taken after 
illumination with the electron beam and the appearance 
of the characteristic diffraction peaks demonstrates that 
the partial measurement of a subset of atoms does not 
destroy the coherence of the remaining system. Further- 
more, it is an example for a complementary measurement 
in position and momentum space on a single many-body 
quantum system. For a quantitative analysis we compare 
the integrated linescans with the Bloch wave function 
that describes the ground state of non-interacting atoms 
in the lattice potential (Fig. 5a, b). The periodic structure 
and the shape of the individual on-site wave function are 
well reproduced for both data sets. Together with the 
observed interference pattern both, the density distribu- 
tion and the quantum-mechanical phase are determined 
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FIG. 5: Ground state of a Bose-Einstein condensate in a one- 
dimensional optical lattice. The graphs show integrated line 
scans for a lattice depth of s = 10 (a) and s — 20 (b), where 
s measures the lattice depth in units of the recoil energy. 
The blue columns are the experimental data which are com- 
pared to a theoretical model (black line) which is based on 
the ground state in the lattice potential (red dashed line) 
convolved with a Gaussian electron beam profile with 95 nm 
FWHM diameter. An absorption image after 15 ms time of 
flight (inset in a) reveals that the phase coherence of the con- 
densate is preserved after exposure to the electron beam. 



jo = i/(27rpo) is the current density in the beam centre, 
/ is the beam current, po is the a- width of the beam, e is 
the electron charge, and <r tot = 9 x 10 -17 cm 2 is the total 
electron scattering cross section for rubidium at 6 keV 
electron energy. For typical beam parameters (I — 23 
nA, FWHM = 140 nm, corresponding to po — 85 nm) 
we obtain r « 17/iS. If the pixel dwell time At is much 
smaller than r, the probability for a scattering event 
(ionization, elastic or inelastic scattering) is given by 
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At j 

— = -<T to t >< A*. 

r e 



If the atom is described by a wave function cf)(x, y, z) 
and if we assume that the beam is much smaller than 
the extension of the wave function, the probability of a 
scattering event at the position y} is given by 

w(x,y) = ^cr tot x At x J dz\cj)(x,y,z)\ 2 . (2) 

Multiplying equation (2) with the ion production effi- 
ciency (Tion/otot, the detector efficiency Tfaet, an d the to- 
tal number of atoms gives equation (1) of the main text. 

Bimodal distribution: For a given number of 
condensed atoms No we numerically solve the Gross- 
Pitaevskii equation 



and thus, the Bloch wave function is fully characterized. 
Eventually, we conclude from the good agreement that 
our imaging technique achieves a spatial resolution of 
better than 150 nm (see Methods section). 

The combination of high spatial resolution and single 
atom sensitivity will open up new possibilities for the 
in situ study of spatial and temporal correlations in 
trapped quantum gases [5j El [7]. Previous experimental 
work on correlations in trapped [27] and expanding 
[4j [19j [28] gases has already demonstrated the high 
potential of such measurements. In optical lattice sys- 
tems, the technique cannot only be used as a powerful 
characterization tool. Removing atoms from specific 
lattice sites, it also allows for the preparation of tailored 
quantum systems. Taking advantage of the electron 
beam's magnetic field even a coherent manipulation of 
single atoms could be feasible. 



Methods 

Detection probability: An atom that is located 
in the centre of a Gaussian electron beam with a 
radial current density of j(p) = jo exp(— p 2 /2p 2 ) ) 1 has a 
lifetime against electron impact of r = e/(jocr to t)- Here, 
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V + Kxt(r"H^(r)| 2 



ip(r) pijj(f) (3) 



using an imaginary time propagation algorithm. The ex- 
ternal potential is cylindrically symmetric and has the 
form V ext (p,z) = \ m {u 2 p 2 + uj 2 z 2 ), where z denotes 
the axial direction of the condensate, u a = 2ir x 12 Hz 
(ujp = 2tt x 170 Hz) is the axial (radial) oscillation fre- 
quency of the dipole trap, p is the chemical potential, 
g = 4:7rh 2 a/m is the coupling constant, and m is the ru- 
bidium mass. For the s-wave scattering length we use 
a value of a = 101 ao, with ao being the Bohr radius. 
In our experiment we produce a spinor condensate in 
the \F = 1) ground state of rubidium. For the model 
presented here we neglect the spinor nature because the 
difference in the scattering lengths for the F — and 
F = 2 scattering channels are only 1% [29]. The con- 
densate wave function is normalized to the total number 
of condensed atoms, Nq = J d 3 x |^(r)| 2 . The numeri- 
cal solution of equation (3) is used to model an effective 
potential for the thermal component 

Kff(r1 = V; xt (r) + 2 5 |^(r)| 2 . 

The density distribution of the thermal component is 
then given by 



n t h(r) 



\ -3 



03/2(2) 



(4) 



5 



with a modified fugacity 

z = exp[-(V e v(r)-fi)/k B T}. 

Here, A t h = \J 2ttTi 2 /mkBT is the thermal de-Broglie 
wavelength, /c# is the Boltzmann constant and T is the 
temperature. The number of atoms in the thermal com- 
ponent is given by 7V t h = / d 3 x n t h(j) and the total num- 
ber of atoms is N = No + 7V t h • 

Spatial resolution: The size of the electron beam 
can be determined independently by scanning the beam 
across a sharp edge of a movable test target which is im- 
plemented in the vacuum chamber. We define the resolu- 
tion as the distance between two neighbouring point-like 
scatterers where the signal intensity in between drops 
to 75%. This definition is the analogue to the Rayleigh 
criterion in optics and for our system translates into a 
resolution of d = 1.18 x FWHM, assuming a Gaussian 
beam profile. The electron beam used for the measure- 
ment in Fig. 4 and 5 of the main text has a diameter of 
95 nm FWHM, corresponding to a resolution of 115 nm. 
The good agreement between the experimental line scan 
and the theoretical model in Fig. 5 proves that a similar 
resolution is achieved for the electron microscope images 
of ultracold atoms. 
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